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The anomaly (cusp) in the mass spectrum of 2^" from Tr^Tr^Tr" decay has been observed in the NA48/2 

experiment. Using the recently developed interpretation of this effect in terms of Chiral Perturbative Theory, the 
pion scattering parameters are measured. The preliminary result for the pion scattering length difference, based 
on the part of statistics, is (ao — a2)m^+ = 0.281 ± 0.007{stat) ± 0.014{syst) ± 0.014:{extern). 



1. Beam and detectors 

The NA48/2 experiment at the CERN SPS is 
searching for the charge asymmetry in the decays 
7r±7r+7r- and tt'^tt^tt". 

The experiment uses two simultaneous oppo- 
sitely charged beams with a central momentum of 
60 GeV/c and a momentum band ±3.8% , prop- 
agating along the same beam line. The decay 
volume is a 114 m long vacuum tank with a di- 
ameter of 1.92 m for the first 66 m, and 2.4 m for 
the rest. The detected particles are dominated by 
the decay products. 

Charged particles are measured by the mag- 
netic spectrometer, consisting of four drift cham- 
bers PP and a dipole magnet, located between 
the second and third chambers. Charged parti- 
cles are magnetically deflected in the horizontal 
plane by an angle corresponding to momentum 
kick of 120 MeV/c. The momentum resolution of 
the spectrometer is a{p)/p = 1.02% © 0.044%p {p 
in GeV/c). 

The magnetic spectrometer is followed by a 
scintillator hodoscope consisting of two planes 
segmented into horizontal and vertical strips and 
arranged in four quadrants. 

A liquid krypton calorimeter |2] is used to re- 
construct tt" — + 77 decays. It is segmented 
transversally into 2x2 cm cells by a system of Cu- 
Be electrodes. The calorimeter is 27 Xo thick and 
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absorbs the electromagnetic showers almost com- 
pletely. Its energy resolution for gamma-quanta 
is a{E)/E = 0.032//B ® 0.09/^; © 0.0042, E in 
GeV. The space resolution for single electromag- 
netic showers can be parametrized as — Cy = 
0.42/-\/E'®0.06 cm for each transverse coordinate 
x,y. 

A neutral hodoscope consisting of a plane of 
scintillating fibres is installed in the liquid kryp- 
ton calorimeter at a depth Q.BXq. It is divided 
into four quadrants, each consisting of eight bun- 
dles of vertical fibres optically connected to pho- 
tomultiplier tubes. 

2. Trigger and events selection 

Data have been collected in 2003 and 2004, pro- 
viding samples of about 10^ K"^ tt'^tt^tt" de- 
cays. Here we report the analysis results of a par- 
tial sample « 2.3 x 10'' decays recorded in 2003. 

7r^7r°7r° decays are selected by a two- 
level trigger. The first level requires a signal in at 
least one quadrant of the scintillator hodoscope 
in coincidence with the presence of energy de- 
positions in calorimeter consistent with at least 
two photons. At the second level, a fast on-line 
processor reconstructs the momentum of charged 
tracks and calculates the missing mass, assuming 
the decay ^ + X of a, 60 GeV/c kaon fly- 
ing along the beam axis. The requirement that 
the mass of missed system X is not consistent 
with the 7r° rejects the main tt^tt^ back- 

ground. 

Events with at least one charged particle hav- 
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ing a momentum above 5 GeV/c, and at least 
four energy deposit clusters in liquid krypton 
calorimeter, each above the energy threshold of 
3 GeV and close in time to the track, are se- 
lected for the further analysis. The distance be- 
tween any two photons in calorimeter is required 
to be larger than 10 cm, and the distance be- 
tween each photon and the impact point of any 
track on calorimeter must exceed 15 cm. Fidu- 
cial cuts on the distance of each photon from the 
calorimeter edges and centre are applied to avoid 
the shower energy leakage. The distance between 
the charged particle track and the beam axis at 
the first drift chamber is required to be larger 
than 12 cm. 

For each possible pair of photons we assume 
that it originates from tt'^ — > 77 decay, and cal- 
culate the distance between the decay vertex and 
the front plane of calorimeter: 



Uik — (,lj 

where Ei, Ek are the energies of the two pho- 
tons and X, y are their impact point coordinates 
on calorimeter. 

The two photon pairs with the smallest Dik dif- 
ference are selected as the best combination con- 
sistent with two tt" mesons from — > Tr^Tr^Tr" 
decay, and the arithmetic average of two D val- 
ues is used as the decay vertex position. The final 
event selection requires that the Tr^Tr'^Tr" invari- 
ant mass differs from the PDG charged kaon mass 
by at most ±6 MeV. 
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3. Cusp effect 

The invariant mass of tt^tt^ subsystem (Afoo) 
in Tr^vr^TT^ decays has been studied to 

search for possible anomaly at the threshold re- 
gion: Moo = 2 • m^±. Figure nia) shows the Moo 
distribution, without any acceptance correction, 
for the selected events. 

Fig^b) shows the same distribution in the 
region close to the threshold Mqq — 4m^± = 
0.0779(Mey/c2)2. The slope change at the 
threshold is clearly visible. Such an anomaly has 
not been observed in previous experiments. 



Figure 1. The M^(7r"7r") distribution for the 
subsystem of the reconstructed tt^tt^tt^ 
decays: (a) ~ in the full kinematic region; (b) 
- in the region around the threshold; the value 
corresponding to 4:Tri^± is indicated by vertical 
line. 



4. Interpretation in terms of CHPT 

The sudden change of slope observed in this 
plot suggests the presence of a threshold "cusp" 
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effect from tfie decay — > tt^tt+tt" contribut- 
ing to the TT^TT^TT^ amplitude through 
the charge exchange reaction tt+tt^ tt^tt". 
The presence of a cusp at Mqq = (2 • m^±)^ in 
7r°7r° elastic scattering due to the effect of virtual 
TT~^Tr~ loops has been discussed first by Meifiner 
et al. [3*. For the case of tt^tt^tt^ decay 
Cabibbo 4^ has proposed a simple re-scattering 
model describing the decay amplitude as the sum 
of two terms representing "unperturbed ampli- 
tude" and contribution from the tt^tt+tt^ 
decay amplitude through tt+tt^ —f tt'^tt*' charge 
exchange. The latter contribution is proportional 
to Oa;, the S-wave tt'^tt' charge exchange scatter- 
ing length. In the limit of exact isospin symme- 
try ttx = (ao — a2)/3, where oq and 02 are the 
TTTT scattering lengths in the 1=0 and 1=2 states, 
respectively. 

Recently Cabibbo and Isidori fjl have proposed 
a more complete model, that takes into account 
all rescattering processes in two-loop approxima- 
tion. In the limit of exact isospin symmetry 
the amplitude depends on five S-wave scattering 
lengths, that can be expressed as linear combina- 
tions of oo and 02 . This model have been used in 
order to extract qq and 02 from the fit to the data. 
The decay amplitude depends not on S3 = Mqq 
only, but also on the two tt^tt^ masses squared 
si,S2. We replace them with the average values 
to avoid multidimensional integration in the fit 
procedure: 



Where 3so = ^Si — rn\ + 2m1a + . 

In order to extract the decay amplitude param- 
eters, the following function F{Mqq) of the mea- 
sured 277° mass was fitted to the experimental 
Mqq spectrum: 

F^N J fiaQ,a2,go,h+,m^MM^Q,m^)dm^ (3) 

Here f is the theoretical distribution from 5. of 
kaon decays over the true value of two-pion invari- 
ant mass m^. $ is the function, describing the ex- 
perimental setup resolution and acceptance, that 
is obtained by Monte Carlo simulation. 



The fit by this model in the interval 0.074 < 
M^o^o < 0.097(Gel^/c2)2 allows to extract four 
parameters: (00 — 02)™^+, a2m^+, qq (the slope), 
(a quadratic term in the Dalitz distribution 
over u variable) and one parameter for normal- 
ization. It leads to = 154.8 for 146 degrees 
of freedom. The better fit (x^ = 149.1/145d./.) 
is obtained by adding to the model a term de- 
scribing the expected formation of a tt+tt^ atom 
("pionium") decaying into tt^tt" at Mqq — 2m^+. 
The fit value for the rate of — > tt^ + pionium 
decays is (1.61 ± 0.66) • 10^^, which is in good 
agreement with the predicted value ^ 0.8 • 10^*^ 
in]. Sensitivity of the measured (oq — a2)mjr+ to 
the pionium production branching ratio is rather 
small: 0.1a(BR)/BR. 

The model used do not include radiative cor- 
rections, which are particularly important near 
Moo = 2mTr+, and contribute to the formation of 
the pionium. For this reason a group of seven con- 
secutive bins with the width of 0.00015(061^/ c2)^ 
centered at Moo — '^1^1^^+ have been excluded from 
the fit in one of the versions of analysis. In such 
a case the measured (ao — a2)m^+ value becomes 
lower by 0.008, that is taken into account as the 
contribution to the systematical error. 

There were made three independent analyses 
which use different event selection criteria and 
different Monte Carlo simulations to take into ac- 
count acceptance and resolution effects. Two of 
them use a fast Monte Carlo code, tuned to the 
experimental resolutions and taking into account 
the geometry and kinematics. For the third anal- 
ysis the full scale GEANT - based [H| simulation 
of the detector is implemented. In all the analyses 
the extensive checks for the accordance between 
the simulation results and data have been made. 
The difference between two simulations leads to 
slightly different second derivative of the accep- 
tance as a function of Moo, that doesn't affect the 
parameter (ao — a2)ni^+. 

Study of systematics shows that no corrections 
are needed. However, the corresponding uncer- 
tainties were taken into account. For measured 
(ao — a2)m7r+ parameter the major ones are shown 
in the tabled 
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Table 1 

Contributions to the systematical error. 

Excluding of pionium region from the fit 0.008 

Cuts variation 0.004 

Z - dependence 0.009 

MC-related difference 0.006 

TOTAL (adding in quadrature) 0.0f4 



5. Preliminary result and discussion 

A critical parameter of the models ^ and [S] is 
the ratio R = ^4+^ /^+oo between weak ampli- 
tudes of TT^TT+TT^ and tt^tt^tt^ de- 
cays. The corresponding "external" uncertainty 
should be complemented by an additional theo- 
retical error of ±5% as the result of neglecting 
higher-order terms and radiative corrections. 

Taking into account these uncertainties the fol- 
lowing preliminary result is obtained: 

(ao - a2)m^+ 0.281 ± 0.007{stat) 
±0.014(syst) ± 0.014(ea:tern) 

This value is compatible withing the errors with 
the results obtained in the E865 BNL j9| and the 
Dirac HU' experiments, as well as with the CHPT 
calculation result jH], ^21 ((o-o ~ a2)™7r+ — 
0.265 ± 0.004). It is also in good agreement with 
the result of fit to low energy experimantal data 
using dispersion relations J3| ((ao — a2)m^+ = 
0.278 ±0.016). 

In the forthcoming stages of analysis we plan 
to decrease the systematical uncertainty and to 
take into account the isospin symmetry break- 
ing. At tree level, omitting the one-photon ex- 
change diagrams, isospin symmetry breaking con- 
tributions to the elastic tttt scattering amplitude 
can be expressed as a function of one parameter 
e — (m7r+ — mJo)/mJ+ ~ 0.065 In par- 

ticular, the amplitudes at the threshold of 271^ 
mass are corrected by a factor of (1 — e) for the 
processes tt+tt"'" Tr~^TT~^, tt+tt^ tt+tt" and 
tt^tt" tt^tt'^; by (1 + e) for tt+tt^ tt+tt^; and 
by (1 -|- e/3) for tt+tt^ — > tt'^tt'^ . Implementing of 
such a correction will lead to the corresponding 
decrease of measured (oq — a2)m^+. 



REFERENCES 

1. D.Bederede et al., Nucl. Instr. and Methods A 
367, 88 (1995). I.Augustin et al., Nucl. Instr. 
and Methods A 403, 472 (1998). 

2. G.D.Barr et al., Nucl. Instr. And Methods A 
370, 413 (1993). 

3. U.-G. Meifiner, G. Miiller and S. Steininger, 
Phys.Lett. B406 (1997) 154: 

U.-G. MeiBner, Nucl.Phys. A629(1998)72. 

4. N. Cabibbo, Phys. Rev. Lett. 93 (2004) 
121801 

5. N. Cabibbo and G. Isidori, JHEP 0503 
(2005) 021 

6. Z.K. Silagadze, JETP Letters 60(1994)689. 

7. G. Colangelo, J. Gasser, H. Leutwyler, 
Phys.Rev.Letters 86(2001)5008. 

8. GEANT: Detector Description and Simula- 
tion Tool, CERN Programming Library Long 
Writcup W5013, version 3.21. 

9. S. Pislak et al. Phys.Rev. D67 (2003) 072004. 

10. B. Adeva et al., hep-ex/0504044;2005). 

11. G. Colangelo, J. Gasser, H. Leutwyler, Phys. 
Lett. B 488 (2000) 261. 

12. G. Colangelo, J. Gasser, H. Leutwyler, 
Nucl.Phys. B 603 (2001) 125. 

13. J.R. Pelatez, F.J.Yndurain Phys. Rev. D71 
(2005) 074016. 

14. K.Maltman and C.E.Wolfe, Phys. Lett. B393 

(1997) 19; Phys. Lett. B424 (1998) 413; 
M.Knecht and R.Urech, Nucl. Phys. B519 

(1998) 329; G. Colangelo, J. Gasser, B.Kubis 
and A.Rusetsky, Comments on isospin break- 
ing in TTTT scattering (unpublished). 



